Stationary solution of a weak-driven open Jaynes-Cummings system of a degenerate 
two-level atom coupled to an arbitrary-polarized cavity field 



On 
O 
O 
(N 

> 

O 

00 



a,: 



> 

oo 
in 
m 

o 

ON 

o 



Sungsam Kang, Youngwoon Choi, Sooin Lim, Wookrae Kim, Jai-Hyung Lee, and Kyungwon 
Department of Physics and Astronomy, Seoul National University, Seoul, 151-742, Korea 

(Dated: October 8, 2009) 

Analytical solution for the stationary density matrix is derived, by using the Morris-Shore transfor- 
mation, for an open Jaynes-Cummings system of a two-level atom with Zeeman sublevel degeneracy 
coupled to an arbitrary-polarized cavity mode. In the limit of weak excitation with the number of 
quantum in the system not exceeding one, we have obtained the stationary solution of the master 
equation up to the first order of the driving field intensity. We have also derived the analytic ex- 
pressions for the excitation spectra of atomic spontaneous emission and cavity transmission. Our 
results show that the system can be regarded as a non-degenerate two-level system with a single 
effective coupling constant which depends only on the elliptic angle of the driving field as long as the 
atom-cavity coupling is not too strong. A precise condition for this approximation is derived. This 
work provides a theoretical ground for experimentally realizing a Jaynes-Cummings system with a 
coupling constant continuously varied for various cavity quantum electrodynamics studies. 

PACS numbers: 42.50.Pq, 37.30.-H, 32.80.Xx 



I. INTRODUCTION 

In cavity quantum electrodynamics (QED), the 
Jaynes-Cummings model (JCM) [Ij is one of the key el- 
ements, describing the interaction between a two-state 
atom and a radiation mode of a cavity. It is not too much 
to say that the cavity QED has been developed within 
the theoretical framework of the JCM. In studying the 
JCM, one is often faced with the difficulty in dealing with 
multi-state atoms in actual problems Q; the JCM itself 
is for a two-state single atom coupled to radiation modes. 
Researchers have proposed extended models for specific 
types of multi-energy-level structure such as three-level 
atoms [1, H[, ladder- type four-level atoms @, two-level 
atoms with Zeeman degeneracy etc. 

Among those multi-state atoms, a two-level atom with 
Zeeman sublevel degeneracy, or the so-called a degener- 
ate two-level system, is of particular interest owing to its 
possible application to quantum state en gine ering [H, 
and quantum information processing [Tol. Illj. Unfortu- 
nately, finding an analytic closed-form description of the 
interaction of the degenerate two-level system with a res- 
onant cavity field is a very complicated matter because 
of the coherence between the Zeeman sublevels. That is 
why this system has usually been studied by numerical 
methods Q except for some limited cases with specific 
polarization d, [TTJ] or specific angular momentum Q ■ 

Nonetheless, as speculated in the review article of Ref. 
2] , it might be possible to find a stationary solution for 
the interaction of a multi-level atom with an arbitrary- 
polarized cavity field in an invariant and analytic form 
with the aid of a unitary transformation introduced by 
Morris and Shore [l^ . For the case of a degenerate two- 
level atom interacting with an arbitrarily polarized light 
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field in a free space, the problem in the semi-classical 
limit, where the atom is quantum mechanical but the 
field is regarded classical, comes down to finding a steady- 
state solution for the generalized optical Bloch equa- 
tion(GOBE) [l3[. Researchers paid a particular attention 
to this problem in order to describe polarization gradient 
cooling [14j | and coherent population trapping with ellip- 
tical dark states [HI], and a full theoretical analysis was 
made [ljl [l7| based on the Morris-Shore transformation 

In this paper, we have extended the previous works 
based on GOBE to an open JCM in a full-quantized man- 
ner. We consider the Zeeman sublevels degenerate in the 
absence of external magnetic fields. The cavity is weakly 
driven by an arbitrary-polarized field while the cavity 
mode is strongly coupled to the F «-> F + 1 closed tran- 
sition of atomic hyperfine structure. We have found an 
analytic expression for the stationary density matrix by 
solving the master equation for the dominant order of the 
driving field intensity. We have also derived the analytic 
expressions for the excitation spectra of atomic spon- 
taneous emission and cavity transmission, respectively. 
Our results show that the system can be regarded as a 
non-degenerate two-level system, unless the atom-cavity 
coupling is not too strong, with a single effective cou- 
pling constant which depends only on the elliptic angle 
of the driving field. We present a precise condition for 
this approximation to be valid. By using our results one 
should be able to realize a JCM with a coupling constant 
continuously varied for various cavity-QED experiments. 

This paper is organized as follows. We formulate our 
system of interest with the master equation in Sec. Inland 
then present a stationary-state solution by introducing 
natural basis in Sec. IIVI The excitation spectra of the 
cavity transmission and the atomic spontaneous emis- 
sion and the resulting non-degenerate two-level system 
are discussed in Sec. |V1 followed by a conclusion in Sec. 

ivm 
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FIG. 1: (a) Left: conventional coordinate frame for elliptic po- 
larization. Right: natural coordinate frame. Polarization vec- 
tor is a superposition of circular and linear polarization, (b) 
Degenerate two-level atom interacting with a common cav- 
ity field, (c) Collective dressed-state representation with the 
number of quantum limited up to one. 

II. THEORETICAL MODEL 

Arbitrary polarization in three dimensions can be ex- 
pressed as 

+1 

e= e i^i W 

q=-l 

where e g 's are unit vectors. Note that in Cartesian co- 
ordinates, these unit vectors are chosen as eo = e Zl 
e± = t(&x ± ie y )/^/2. It is always possible to remove 
one component of the polarization vector through an ap- 
propriate choice of e z axis. For example, if we choose 
e z axis to be perpendicular to the plane of polarization, 
the e z component vanishes, and the polarization becomes 
elliptical in general as shown in the left part of Fig. [Ha). 

Throughout this paper, we choose a natural coordinate 
frame as introduced in Ref. The unit vectors e£ of 

this coordinate frame are shown in the right part of Fig. 
[IJa). In the natural coordinate frame, one of the circular 
component vanishes, and thus the polarization vector is 
a superposition of a circular and a linear polarization 
component: 

e = egV cos 2e — e' +1 v2 sin e, (2) 

where e is the elliptic angle defined by Fig.QJa), satisfying 
— 7r/4 < e < +7r/4. In this frame, linear polarization 
can also be described naturally (e = 0), and moreover, 
the analytic expressions of the matrix elements for the 
stationary density matrix become greatly simplified as 
to be shown below. 

Our system of interest is shown in Fig. [IJb) . A degen- 
erate two-level atom is coupled to a cavity with a repre- 
sentative coupling constant go while the cavity is weakly 



driven by a classical field of arbitrary polarization e and 
an amplitude £. The decay rates of the cavity and the 
atom are denoted by k and 7 (both half widths), respec- 
tively. We assume that the frequency of the driving field 
(lul) and that of the cavity (u>c) are both near resonant 
to the atomic transition (oja), (F g = F) «-> (F e = F + 1). 
The coupling constant 2go is the vacuum Rabi frequency 
for the cycling transition, (mj? = F) <-> (m' F = F + 1). 

Under the rotating wave approximation, the Hamilto- 
nian for this system can be described by a multi-level 
expansion of the Jaynes-Cummings Hamiltonian Q . In 
the rotating frame of ujl {i.e., in the interaction picture), 
the Hamiltonian is written as 

+1 

H/h = A c a^a + A A £ D\D q 

q=-l 

+ g (Va^ + aV^)+£(a + a^) , (3) 

where Aa,c — {^a,c~^>l) are the detunings of the atom 
and the cavity with respect to the driving field, respec- 
tively. The classical-driving-field amplitude £ is scaled 
in such a way that £/k represents a dimensionless in- 
jected photon flux. Operator D q is defined as an atomic 
lowering operator for polarization 

D *= E C F g ili q \a,rn F ){e,m' F \ , (4) 
with the Clebsch-Gordan coefficients 

c^ m ^(~iy-^ M V2JTT^ 1 _ J M ) . (5) 

Note that D q satisfies the relation F>\D q = 
J2 m ' l e ' m F)( e > m F\i which implies the conservation of 
the total population for a closed transition. Operator V 
is an atomic lowering operator for e polarization and can 
be written as 

+1 

V = ( E D «K) ■ e (6) 

g=-l 

Since mp = —F,—F+1,...,F, and m' F = — F — 
1, — F, . . . , F + 1, both operators D q and V can be repre- 
sented by (2F+1) x (2F+3) matrices. The time evolution 
of the system is then described by the master equation, 

p = C[p] = —[H, p] + K(2apa^ — a^ap — pa^a) 
in 

+1 

+ 7E ( 2D «P D t - D l D q P - pD\D q ), (7) 

q=-l 

and the stationary density matrix is obtained by solving 
C[p ss ] =0. 
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III. WEAK EXCITATION LIMIT IN 
COLLECTIVE DRESSED-STATE BASIS 

In general, an analytic expression for the stationary so- 
lution of Eq. (J7JI cannot be obtained because the field op- 
erators a and a' have an infinite number of bases. How- 
ever, in the weak excitation limit, one can expect that 
the infinite number of field bases can be truncated to a 
finite number of low-quantum bases and that the steady- 
state solution can be expressed approximately in terms 
of a few low orders of driving- field amplitude £ [19|, l20| • 
Under this assumption, the system can be described by 
the interaction between two manifolds of (2F + 1) collec- 
tive dressed states |a) co i and (2F + 3) collective dressed 
states 1 6) col with the number of quantum for both to be 
one. These manifolds are connected to a manifold of 
(2F + 1) collective ground dressed states | 0) co i with no 
quantum through two decay channels as shown in Fig. 
[]Jc). These collective state vectors are defined as 

|0)col = Em*. £m P lf= ™F>atom <g> 1 0) field (8a) 
|a) co l = E mF ^m P \9,m F ) atom ® |1) field (8b) 
|&)col = J2rn' F d' F l e ' TO f)atom ® |0) fie ld , (8c) 

where {£j°' a } are complex amplitudes. There are as 
many as (2F g ^ e + 1) independent sets of {^j °' a,f, - ) }. In the 
absence of the driving field {£ = 0), only can have 
non-zero values in the steady state. Therefore, we expect 
that the dominant terms of b } are at least of the first 
order of £ while are of the zeroth order of £ . 

With these bases the Hamiltonian and the density ma- 
trix can respectively be expressed by a (6F+5) x (6F+5) 
matrix, which can be divided into 9 sub-matrices, as 



H 



g V AchF+i 
£I 2 f+i 



Pbb Pba PbO 
Pab Paa PaO 
POb POa POO 





£I2F+1 




(9) 



(10) 



where /„ denotes an n x n unit matrix and pij = 
I i) col 01 col for {i, j} = {0,a, b}. By substituting Eqs. © 
and (|10[) into Eq. ([7]), we obtain a set of first-order dif- 
ferential equations 

p bb = -2-ip b b - igo(V^pab ~ PbaV) (11a) 

Paa = -2np aa -ig (Vpba~ PabV f )~i£(poa- PaO) (lib) 

poo = 2np aa + 27 ^ D q p bb D q + i£(p Qa - p a0 ) (11c) 

Pab = {pba) ] = i{E* A -E C )Pab-igo{Vpbb- PaaV)-i£pob 

(12a) 



Pab = (pfco) f = iE*iPob + igoPoaV - i£p a b (12b) 

POa = (paoV = iEcPQa+igopobV^ +i£{p o- paa), (12c) 

where Ea = A A — «7 and Ec = Ac — in are the complex 
energies of the atom and the cavity, respectively. 

To solve for a stationary solution, we should find 9 sub- 
matrices which make the time derivatives in the left-hand 
side of Eqs. (fTTj) and ([12")) vanish, under the normalization 
condition Tr[p] = 1. By rearranging terms and removing 
the higher order terms of £ than the second order, we 
obtain 

Pab = (Pba) 1 = is*g (Vpbb - PaaV) (13a) 
+ is*g £ 2 pooVX- 1 + 0(£ 4 ), 



Pob = {pbo) f = g £pooVX- 1 + 0(£ a ), 



(13b) 



Poa = (PaoY = -E A £p m Y~ l + 0(£ 3 ), (13c) 

where s = i/ {Ea — Eq). X and Y are the square matrices 
defined as 

X = E* A E* C I 2F+3 - glV^V , 
Y = E* A E* C I 2F+1 - g 2 VV^ . 

One can easily prove X and Y satisfy the following iden- 
tities: 



XV^ = V^Y, VX = YV. 



(14) 



By substituting Eqs. |(T3j! and (O into Eq. (JTTJ, we 
then obtain 

~2 1Pbb - [s*glV\p aa -£ 2 p m Y- x )V 

-s*g 2 V Wp bb + h.c] =0, (15a) 

- 2npaa + [s*g 2 Q (paa ~ £ 2 PooY^VV^ 

~s*glVpb b V J< +iE* A £ 2 Pm Y- 1 + h.c] = 0, (15b) 

2KPaa + 27 ^ D q p bb D\ 
9 

-[iE A £ 2 p QQ Y- 1 +h.c]=Q, (15c) 

where 'h.c' denotes the Hermitian conjugate. The Eq. 
(Tl~5|) is still too complicated to solve. It is noted, however, 
that except for the second term of Eq. (|15cp . all terms 
are connected only by the dipole transition matrices V 
and V ' . Therefore, one may expect that Eq. (|15[) can be 
greatly simplified if the density matrices can be simul- 
taneously diagonalized with both V and V* . Note that 
both V and are not square matrices. The meaning of 
diagonalizing these non-square matrices is to make only 
the diagonal components of these matrices have non zero 
values as specified in Eq. (jTTJ) below [l2T |. 
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T~ 1 

^ 1 ^2F+1 

(a),l) \(a),2F + \) 

FIG. 2: (a) Conventional Zeeman basis and transition link- 
age by polarization e. (b) Natural basis and the correspond- 
ing transition strength A^. Two sublevels of the upper state 
manifold are uncoupled. 



IV. STATIONARY-STATE SOLUTION OF 
MASTER EQUATION IN NATURAL BASIS 

The system of a degenerate two-level atom interacting 
with an arbitrary-polarized classical field can be trans- 
formed to a set of independent non-degenerate two-level 
systems by an unitary transformation, known as the 
Morris-Shore transformation [ijj ■ As shown in Ref. [2l| , 
VV^ and V'V become diagonalized simultaneously un- 
der this transformation and can be written as 

2F+1 

Vyj = E X l(a)Mi)((*)il 

i=l 
2F+3 

V ^ V = E #m\Q>)3)((P)3\- (16) 

3=1 

The states \(a)i) and \(b)j) are the eigenvectors of VV^ 
and V^V, respectively, with X 2 ^ and A^/^ correspond- 
ing eigenvalues, respectively. It can be proven that the 
eigenvalues are real non-negative numbers and that the 
sets {A?( a )} and {A 2 , b -.} of non-zero eigenvalues coincide 
with each other. 

Accordingly, as shown in Ref. [2l| . we can write the 
dipole interaction operators V and in a diagonal form 
in this new basis, called the natural basis [2ll |. as 

V = J2M(a)i){(b)i\, = '£\ i \(b)i){(a)i\. (17) 

i i 

where Xi,( a ) — ^i,(b) — Aj for i = 1,2,..., 2F + 1 and 
XjJb) — for j = 2F + 2 and 2F + 3, which implies that 
no transition is allowed for \(b)2F + 2) and \{b)2F + 3) 
sublevels [see Figl^b)]. 



In the conventional Zeeman basis, all degenerate sub- 
levels are linked together by dipole transitions [see 
FigEfa)]. In the natural basis, however, one ground-state 
sublevel \(a)i) is coupled to only one excited-state sub- 
level \(b)i) with a transition strength A^ [see Figl^b)], like 
a transition by linear polarization in the Zeeman basis. 
Indeed, for linear or circular polarization, the Zeeman 
basis coincides with the natural basis. 

It is shown in Ref. |2l| that we can uniquely determine 
two Hcrmitian matrices A and B which satisfy 

AV = VB, A = J2 q D q BD\ (18) 
[A, VV^} = 0, [B, VW] = 0. (19) 

The stationary-state solution of the density matrix can 
then be expressed with these matrices. From Eqs. (|17[) - 
([12]) we can show that these two matrices are also diag- 
onal in the natural basis with the same set of non-zero 
eigenvalues just like V and V'. Therefore, in the natural 
basis these matrices are expressed as 

2F+3 2F+1 

B= E vj\(P)j)((m> -4= E "il(o)i><(<0*1, ( 20 ) 

3=1 i=l 

where fi's are eigenvalues of A for i = 1,2, ... , 2F + 1 
and v 2 f+2 = v 2 f+3 = 0. 

From Eqs. @, (H]) and ([5]), it can be seen that ma- 
trices V and V' and thus A and B depend only on the 
polarization (elliptic angle e) and the atomic level struc- 
ture (angular momentum F). They do not depend on 
the other parameters such as detunings, intensity, cou- 
pling constant, etc. Moreover, although we are dealing 
with F <-> F + 1 transition, the matrices A and B can 
be found for arbitrary transitions. The details of finding 
these two matrices are well documented in Ref. [2l| and 
references therein, and therefore, in this paper we only 
give the explicit expressions for matrix elements of A and 
B in Appendix lAl 

Substituting Eqs. (JTHJ) and (JT5J) into Eq. (T5J and after 
some straightforward algebra with an assumption that 
Pbb is proportional to B, we finally obtain the diagonal 
elements of the density matrix in closed forms as 

Pbb = V £ 2 g 2 B + 0(£ 4 ) (21a) 
p aa = r 1 £ 2 \E A \ 2 (VV^)- 1 A + 0(£ 4 ) (21b) 

P oo = v(YY^(VV^)- 1 A+0(£ 2 ), (21c) 

where rj is a normalization constant. The results show 
that pbb is indeed proportional to B, which is self- 
consistent with the above assumption, and that all the 
sub-density matrices pij are diagonal in the natural ba- 
sis. Furthermore, it is possible to prove that this solution 
is unique in the same way as described in Ref. [l6l |. 

As mentioned previously, the dominant order of £ is 
£° for poo, and £ 2 for p b b and p aa . The normalization 
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constant r\, determined by the condition Tr[p] — 1, is 
given by 



\E A E c \ia - 2Re[E A E c ]a l9 l + a 2 gt 



Q(£ 2 ), 



with coefficients 



(22) 



a a = Ur[(VVt)-U], a 2 = Th:[VVU], 

ai = Tr[^] = Tr[S]. (23) 

The analytic expressions for these coefficients are given 
in Appendix [Bl 



V. EXCITATION SPECTRA AND 
EQUIVALENCE TO A NON-DEGENERATE 
TWO-LEVEL SYSTEM 

In the previous section, we have derived the ana- 
lytic expression for the stationary density matrix for the 
system of a degenerate two-level atom coupled to an 
elliptical-polarized quantized cavity field. We can now 
obtain the information on the emission from the system 
from this stationary density matrix. There exist two out- 
put channels in our system as shown in Fig. [1] The cav- 
ity transmission (T cav ) is characterized by k while the 
atomic spontaneous emission (T sp ) is characterized by 7. 
These two types of emission can be described by using 
the results from Eqs. l[2"T ]) -(|2"3" ]) as 



T cav = Tr [a f ap] = Tr[p aa ] 



£ 2 \E A \< 



\E A E C \ 2 - 2Re[E A Ec}(a 1 /a )g 2 + {a 2 /a )g* 



T — 



+0(£ 4 ), 

Tr[J2DlD g p]=Tr[p bb ] 
1 

£ 2 (a 1 /a )g 2 



(24) 



\E A E C \ 2 - 2Re[E A E c ](a 1 /ao)g 2 + (a 2 /a Q )g l 



-0(£ 4 



(25) 



They are the expressions for the excitation spectra for 
the cavity and the atomic emission channels, respectively, 
i.e., each describes the normalized emitted power from 
the corresponding output channel as a function of the 
driving laser frequency. 

When the field is circularly polarized (e = ±7r/4), the 
two spectra coincide with those of the non-degenerate 
two-level case [HI!!, by using Eq. ([B6j) . 



T cav (e = ±~)=£' 



T sp (e = ±-)=£ : 



E A 



E A E C 



,9o 



9l 



E A E C 



9l 



+ 0(£ 4 ), (26) 



+ 0(£ 4 ). (27) 



Moreover, the structures of Eqs. (|24[) and (|25|) look sim- 
ilar to those of the non-degenerate two-level case. By 
letting 



.9 



(28) 



the denominators of Eqs. (12"4")) and ([23)1 can be written as 

\E A E C - g' 2 \ 2 + 5(e)g' 4 , (29) 

where (5(e) = aoa 2 /a 2 — 1, which is much smaller than 
unity as to be shown below. Only when the last term 
of the above expression is negligible, the excitation spec- 
tra for two output channels become of the same form as 
those of the non-degenerate two-level system with g' in- 
terpreted as an effective coupling constant of the system. 

We can find the precise condition under which we can 
neglect the last term of Eq. ([29]) . The first term of Eq. 
([29]) can be expanded by using the definition of E A ^c as 



\E A E c -g l2 \ 2 



(A C A A - g 12 ) 2 + (A AK ) 2 
+ [(«7 + 5' 2 ) 2 



.9 



(Ac7) 2 
(30) 



The lower bound of the above expression is 



\E A E c -g' 2 \ 2 >4g' 2 K1 , 
|^c- 5 ' 2 | 2 >(.9 ,2 + «7) 2 , 



if k 7 < g 12 , (31) 
if k 7 > g' 2 . (32) 



Both are combined to a single condition regardless of the 
value of K'f/g' 2 as 



\E A E c -g' 2 \'>4g u K j 



(33) 



from which we can deduce that the condition for neglect- 
ing the last term of Eq. (f2"9")) is 



(34) 



g' 2 9o 



Here we use the inequality go > g' from Eq. (|B10j) . 



VI. DISCUSSION 
A. Accuracy of Two-Level Approximation 

In Table [H we present the maximum value of 5(e) for 
various F with relevant atomic species. The variation 
of 5(e) as a function of e is shown in Fig. [3]Ja) for some 
F values. For circular polarization, 5(e) becomes zero, 
consistent with the fact that the system becomes a non- 
degenerate two-level system in the steady state due to 
optical pumping. We also present the effective coupling 
constant g 1 as a function of elliptic angle e in Fig. E^b). 
As mentioned before, when e = ±7r/4, g' coincide with 
go, and g' has its minimum value when e = 0. 

Since 5 is of order of 0.01, the condition in Eq. (|3~4"1) 
can be safely satisfied if An^f/g^ is order of 0.1 or larger. 
This requirement can be met in a weak coupling regime 
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it/8 7t/4 

elliptic angle, e 



it/4 
elliptic angle, e 




FIG. 3: (Color Online) (a) Plot of 6(e) for various F. (b) 
Variation of g' /go as elliptic angle e. (c) Excitation spectrum 
of the cavity transmission for linear polarization (e = 0) when 
F — 3. For a fixed value of 7/(70 = 0.1, the value of n/go 
is varied. Solid line represents the spectrum given by Eq. 
(|24[) while the dotted line shows the equivalent non-degenerate 
two- level system with a coupling constant equal to g' . 



TABLE I: The maximum values of 5(e) for various angular 
momenta F and corresponding atomic species. 



F 


<5max 


relevant atomic species" 


1 


0.043 


l H 


2 


0.027 


7 Li, 23 Na, 39 K, 41 K, 87 Rb 


3 


0.016 


85 Rb 


4 


0.010 


133 Cs 



a D2 line of these atoms have F *-* [F + 1) closed transition. 



as well as in intermediate and strong coupling regimes. 
Some example are given below. If the coupling is too 
strong, however, the condition is not satisfied. 

In Fig. 03c), the excitation spectrum of the cavity 
transmission, given by Eq. of the degenerate two- 

level system when driven by linear polarization (e = 0) 
is compared with that of the non-degenerate two-level 
system whose coupling constant is the same as g' of the 
degenerate two-level system. Both exhibit a double-peak 
structure, also known as the normal mode splitting. For 
the same 8 = 1 and 7/50 = 0.1, the two spectra are 
compared for three different values of k when F = 3, for 
which 6(0) ~ 0.01. When n/g = 0.01 (i.e., too strong 
coupling), the condition in Eq. (f34|) is not satisfied, and 
consequently, we see a noticeable discrepancy between 
the two spectra (red solid line vs. red dotted line). For 
other values of n in Fig. EIc), however, the condition is 



well satisfied, and the two spectra also agree well with 
each other. 



B. Physical Interpretation of Effective Coupling 
Constant g' and Small parameter 5 

The analytic expressions for g' and S(e) can be derived 
explicitly by using Eqs. (|Bip ~ (|B3|) . However, we gain 
more physical insights by expressing g' and 6(e) in the 
natural basis. With substitution of Eqs. (fll))) and (p?0"|) in 
Eq. (|21b[) . we can express p aa in terms of the populations 

in the natural basis. We define ir^ to be the ratio of the 
population in \(a)i) sublevel to the total ground-state 
population p aa : 



a) _ Vil>% 



(a) _ 



(35) 



1 ~E«"iA?' 

From Eqs. ([To]). (|20f and (J23J) , we have 

a = ^2 Vi/^i, Qi = ^ a 2 = ^2 v i^i- ( 36 ) 



Then, the g' , and 6(e) can be rewritten as, 



fJ 




(37) 



6(e) = a a 2 /al - 1 

(E 4 A^ ) )~(E,AN a; 



(AA 2 ) 
A 2 



(38) 



where x = J2i n i x i the population- weighted average of 
x in the natural basis and (Ax) 2 = x 2 — x 2 the variance of 
x. Therefore, the effective coupling constant g' is scaled 
to go by the root-mean-square of the transition strength 
Xi while the mean is given by the population-weighted 
average in the natural basis. Likewise, 6(e) is the ratio 
of the variance to the square of the population-weighted 
average of transition probability X 2 . In the statistics ter- 
minology, 6(e) is the same as the square of the coefficient 
of variation for the transition probability. With this in- 
terpretation of 6(e), we can have an alternative under- 
standing why (5(±7r/4) = for circular polarization: it is 

due to the vanishing dispersion AA 2 since only one ir^ 
is nonzero in this case. 

As shown in Figs.[3{a) and[3Jb), 6(e) decreases as the 
angular momentum F increases while the change of g' is 
negligible. This means that the variance of the transi- 
tion probability decreases as the F increases whereas the 
average of it does not change much. This can be under- 
stood by considering the case of linear polarization, in 
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FIG. 4: Plot of Af and 7 r l (a) for various F values when e = 0. 
As F increases, the difference of Af between three most- 
contributing sublevels (i = —1,0,1) decreases, and accord- 
ingly, the dispersion AA 2 decreases. 



which the dependence on F is pronounced. When e = 0, 
the transition probability Af becomes, 



A 



(2F + 2-i)i 
(F + 1)(2F + 1)' 



(39) 



with i — 1,2,..., 2F + 1. It has a maximum value 
(F + 1)/(2F + 1) for i = F + 1, which corresponds to 
m.F = in the Zeeman basis. The difference between 
the maximum value and the second maximum value is 
1 — l/(F + l) 2 , which converges to unity as F increases. 
Likewise, the difference between the maximum value of 
izf" and the second maximum value of 7rj also con- 
verges to unity as F increases. Therefore, the variance 
of the transition probability decreases as F increases. 
Fig. 0] shows the plot of Af , and ir^ as F is varied. 
As F increases, the difference of Af between three most- 
contributing sublevels (i = —1, 0, 1) decreases. 



C. Possible Application to Cavity-QED 
Experiments 

We have shown that the steady state properties of a de- 
generate two-level Jaynes-Cummings system in the weak 
excitation limit is approximately the same as that of a 
simple non-degenerate two-level system with an effective 
coupling constant. This result can be applied to many 
cavity QED experiments which should consider linear or 
elliptic polarization. For example, if one employs special 
cavities such as microspheres 1221. microtoroids [23| and 
photonic crystal microcavities [24j |. the optical pumping 
to a cycling transition cannot be fulfilled because circu- 
lar polarization is not supported there due to the special 
geometry of supported modes. Therefore, it is necessary 
to consider the multi-sublevels of atom and coherence 
among them, making the analysis very complicated. One 



can still use linear polarization with these cavities and by 
using our result the original system can be transformed 
to a system of a non-degenerate two-level atom with a 
single effective coupling constant. 

Another important point to note is that the effective 
coupling constant can be a new parameter which can be 
continuously controlled. So far, there has been no way to 
control the coupling constant continuously in the cavity 
QED experiments with a single atom. With our result, 
however, the effective coupling constant can be contin- 
uously changed simply by varying the elliptic angle of 
a probe laser. In particular, with appropriate choice 
of go and k, it should be possible to control the effec- 
tive coupling constants from the weak coupling to the 
strong coupling regime, and thus one can study the quasi- 
eigenstates of the atom-cavity composite near the transi- 
tion point from the weak to strong coupling regime. Such 
studies are of considerable interest related to an excep- 
tional point (2f| and its associated branch point singu- 
larity in the parameter space of the system described by 
a non-Hermitian Hamiltonian 12611. 



VII. CONCLUSION 

We have obtained an analytic solution for the station- 
ary density matrix of an open Jaynes-Cummings sys- 
tem composed of a degenerate two-level atom and an 
arbitrary-polarized cavity field. The stationary density 
matrix is obtained up to the first-order of driving-field 
intensity. We employ a natural basis formed by the eigen- 
states of the dipole interaction matrices and show that 
the stationary density matrix is diagonal in the natu- 
ral basis. We have also calculated the excitation spectra 
for the cavity transmission and the atomic spontaneous 
emission, respectively. Unless the atom-cavity coupling 
is too strong, the precise meaning of which we specify, the 
excitation spectra are shown to be the same as those of 
an equivalent non-degenerate two-level system with an 
effective coupling constant g' . We show that the effec- 
tive coupling constant is a function of elliptic angle and 
angular momentum only, not depending on the other pa- 
rameters. 

The present work is the first to report a full-quantum 
mechanical analytic solution for the cavity QED system 
of an degenerate two-level atom with arbitrary polariza- 
tion under weak excitation. Our results provide a use- 
ful theoretical background for better understanding this 
versatile system and a design basis for future quantum 
optics experiments including quantum information pro- 
cessing and quantum state engineering. Furthermore, the 
fact that the effective coupling constant can be continu- 
ously controlled suggests a new venue in the cavity QED 
experiments such as the study of an exceptional point 
and its associated singular properties in the atom-cavity 
composite described by a non-Hermitian Hamiltonian. 
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APPENDIX A: EXPRESSION FOR MATRIX 
ELEMENTS OF A AND B 

The matrix elements of A and B are determined only 
by the elliptic angle e and angular momentum F. For 
F <-> F + 1 transition, the elements of these matrices can 
be calculated by using Eq. (B8) of Ref. [2l| as 



A j 



Min[i,j" 

E 



i 



i (i-ky.(j-k)\ 



' (2F + 1 + i - fc)!(2F + l+j-k)\ 
(2F+l-i + k)\(2F + 1 - j + k)\ 

sine \*+i-2fc 



'F i (F+l) -fc^F j (F+l) -fc ^ ^/^^ 



(Al) 



Si 



Min[i,j] 

E 

k=-F 



(i-k)\(j-k)\ 



\2F + l + i- k)\{2F + l+j-k)\ 
(2F + 1 - ?: + k)\(2F +l-j + k)\ 



/ sin e \ 



i+j-2k 



(A2) 



The dimension of A is (2F + 1) x (2F + 1) while S is 
(2F + 3) x (2F + 3). 



APPENDIX B: ANALYTIC EXPRESSIONS FOR 
COEFFICIENTS a , on, AND a 2 



The coefficients o>q, a\, and ai in Eq. (|23|) depend 
on elliptic angle e and angular momentum F only. For 
arbitrary polarization e, as defined by Eq. and for the 
transition F <-» F + 1 , these coefficients can be expressed 
in explicit form by using Eq. (79) of Ref. [2l[ as 



1 



a = 



F 

— y 

cos 2e f-^ 



C 2 2 i P 2/ (l/cos2e) 



(Bl) 



Q'l = 



2F+1 



(l/cos2e) 



(B2) 



2F+2 

a 2 = cos 2e ^ D?Pi (1/ cos 2e) , (B3) 

i=2F 

where the Pi(x) denotes the Legendre polynomials, and 
the coefficients Ci , and Di are given by Eq. (77) of Ref. 
(H as 



(2Z + 1)(2F-Z)!(2F + J + 1)! 
G - V (2F+1)(4F+1)! (B4) 



A = v/(2F + 3)(4F + 3)Cl° 02F+10 

/ 1 2F + 1 
F F+ 1 F + l 



(B5) 



where the last factor denotes the 6j symbol. From 
1/ cos2e > 1, we note that ao, cn\ and a 2 should have 
positive values. When e = ±7r/4, these coefficients di- 
verge to infinity, but the ratios between them converge 
to unity: 



a 2 



lim — = lim 



(B6) 



From the recursion relation of Legendre polynomials 

(21 + l)xPl(x) =(l + l)F +1 (x) + lPi-x{x), (B7) 
the coefficients ao and a 2 in Eqs. (|Bip - (|B3[) become 
2F 



ao = a\ 



2F+ 1 



-F 



2F-1 



(l/cos2e) 



F-l 



+ ^E C 2^(l/cos2e), (B8) 



1=0 



a 2 = ci'i 



2Fcos2e 
4F+ 1 



F 2 f (1/ cos 2e). 



(B9) 



From the above equations, these coefficients satisfy the 
following inequality 



ao > ai > a2, 



(BIO) 



and the equality holds for e = ±7r/4 as shown in Eq. 
(IB61). 
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